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ABSTRACT 

Context. The Australia Telescope Large Area Survey (ATLAS) aims to image a 7 deg 2 region centred on the European Large Area ISO 
Survey - South 1 (ELAIS-S1) field and the Chandra Deep Field South (CDF-S) at 1.4GHz with high sensitivity (up to a ~ 10 /Jy) 
to study the evolution of star-forming galaxies (SFGs) and Active Galactic Nuclei (AGN) over a wide range of cosmic time. 
Aims. We present here ancillary radio observations at a frequency of 2.3 GHz obtained with the Australia Telescope Compact Array 
(ATCA). The main goal of this is to study the radio spectra of an unprecedented large sample of sources (~ 2000 observed, ~ 600 
detected in both frequencies). 

Methods. With this paper, we provide 2.3 GHz source catalogues for both ATLAS fields, with a detection limit of 300 pJy (equivalent 
to 4.5<r in the ELAIS-S1 field and 4.0cr in the CDF-S). We compute spectral indices between 1.4 GHz and 2.3 GHz using matched- 
resolution images and investigate various properties of our source sample in dependence of their spectral indices. 
Results. We find the entire source sample to have a median spectral index a mei = -0.74, in good agreement with both the canonical 
value of -0.7 for optically thin synchrotron radiation and other spectral index studies conducted by various groups. Regarding the 
radio spectral index as indicator for source type, we find only marginal correlations so that flat or inverted spectrum sources are 
usually powered by AGN and hence conclude that at least for the faint population the spectral index is not a strong discriminator. 
We investigate the z~a relation for our source sample and find no such correlation between spectral index and redshift at all. We do 
find a significant correlation between redshift and radio to near-infrared flux ratio, making this a much stronger tracer of high-z radio 
sources. We also find no evidence for a dependence of the radio-IR correlation on spectral index. 

Key words, catalogues - galaxies: active - galaxies: evolution - radio-continuum: galaxies - surveys 



1. Introduction 

Large-scale radio surveys are the primary instrument in the tool- 
box of contemporary astronomy to investigate demography and 
evolution of active galaxies across cosmic time. Consequently, 
there have been numerous attempts to use this tool, whether 
in a wide, shallow layout (e.g. the Nation al Radio Astronom 
Observatory Very Large Array Sky Survey. ICondon et"aDll99 
or in a deep, but more pencil beam-li ke layout such as for the 
Cosmic Evolution Survey (COSMOS, Schi nnerer et alJl2010h or 
even the Lockman Hole survey bv lOwen & Morr ison (2008|), the 
most sensitive radio observation to date. 

The Australia Telescope Large Ar ea Survey (ATLAS, 
iNorris et al.l2006l; lMiddelber g et alj2008l) attempts to break this 
dichotomy by obtaining deep images (<x ~ 10/iJy in the fi- 
nal data release 3) over an unprecedented large area of sky. 
Distributed over the well-studied Chandra Deep Field South 
(CDF-S) and the European Large Area ISO Survey - South 1 
(ELAIS-S1) field, the total area covered by ATLAS is approxi- 
mately seven square degrees. This allows us to stud y the faint ra- 
dio sky without being affected by cosmic variance dMoster et alJ 



1201 ll) . Furthermore, the wealth of ancillary photometric data 
ranging across the entire elect romagnetic spectrum (e.g. from 
SWIRE, COMBO- 17 or GEMS lLonsdale et alJ2003l:IWolf et al.1 
2003; Rix et al. 2004) as w ell as numerous spec t roscopic cam- 
paigns (ISacchi et all 12009b Popesso eTaT] 120091: Balestra et al. 
2010; Ma o et alJl2012h enables one to conduct detailed studies 
of radio-detected objects. 

In this paper, we focus on radio spectral properties of 
the ATLAS source sample. Th e m ain ATLAS data r elease s 
are DR 1 by INorris et alJ d2006l) and iMiddelberg etail <|2008), 
DR2 by Hales et al. (in prep.), and DR3 by Banfield 
et al. (in prep.). These data consist of 1.4 GHz obser- 
vations as for most contemporary radio surveys, and we 
here supplement these data with observations at 2.3 GHz. 
The main goal of this is to calculate spectral indices a = 
log(Si.4GHz/S2.3GHz)/log(1.4GHz/2.3GHz) for all sources 
detected at both wavelength and so to characterise their radio 
emission mechanism, source type and dependence on other spec- 
tral quantities. Even though many blind surveys of extragalac- 
tic radio sources are available already, only few studies have 
been carried out which characterise the radio spectra of a large 
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sourc e population by m ea ns of their spectra l index. For exam- 
ple, |BondLitaD (120071) and llbar et all (l2009h use supplementary 
6 10 MHz data to calculate spectral indices at frequencies below 
1.4 GHz with the main aim to identify steep spectrum sources 
which are regarded as good candidates for being high-redshift 
objects (see Sect. 14.3b . The most comprehensive study of spec- 
tral indice s between 1 .4 GHz an d a higher frequency was con- 
ducted by Pra ndoni et al.l d2006l) using 5 GHz Very Large Array 
(VLA) data in the Australia Telescope ESO Slice Project survey 
field. With a sensitivity of ~ 70 /iJy over an area of one square 
degree, their source sample consists of approximately 100 ob- 
jects which are detected at both frequencies, so our sample ex- 
ceeds this by a fa ctor of six. 

According to lKomatsu et alj d2011l) . we adopt a flat ACDM 
cosmology with ffo = 70.2kms~ 1 Mpc~' and Qa = 0.725 
throughout this paper. 

2. Observations and imaging 

In the following description, we give the actual numbers for 
all things concerning the ELAIS-S1 field first, followed by the 
numbers for the CDF-S in parentheses. If there is only a sin- 
gle number given, it is approximately equal for both fields. 
Observations at 2.3 GHz for the ATLAS fields were carried out 
between December 2006 and March 2007 on 52 separate days 
where the ATCA was only in its 750B and 750C configurations. 
This resulted in a substantially lower resolution than at 1 .4 GHz 
since the longest baselines used for imaging at 2.3 GHz was 
750 m, or 5.7 kA, resulting in an angular resolution of 33" (57"), 
compared to 30kvl at 1.4 GHz. The ATCA primary beam at a 
frequency of 2.3 GHz has a FWHM of 21 arcmin, therefore 76 
(88) overlapping pointings were observed in a hexagonal lattice 
with an average integration time of 6.17h (4.13h) per pointing 
(compared to approx. 12 h per pointing for the 1.4 GHz image) 
to cover the entire area of the 1 .4 GHz maps. 

For phase calibration, which was done approximately ev- 
ery 45 min, PKS 0022-423 (PKS 0237-233) was observed as sec- 
ondary calibrator. Flux calibration was done for both fields us- 
ing PKS 1943-638 which was observed at least once a day 
for 10 min. As for 1.4 GHz, the ATCA provides two indepen- 
dent frequency bands with a bandwidth of 128 MHz divided 
into 33 channels. During the first two days of observations for 
the ELAIS-S1 field, the two bands were centred at 2.368 GHz 
and 2.560 GHz in an attempt to reduce local Radio Frequency 
Interference (RFI). The fluxes density of the primary calibrator, 
PKS 1934-638, was assumed to be 1 1.589 Jy and 10.932 Jy, re- 
spectively. Because using a higher frequency for the upper band 
did not significantly decrease the RFI, the remaining time was 
observed with the standard band configuration, meaning that 
the upper band was lowered in central frequency to 2.496 GHz. 
There, the primary calibrator was assigned to have a flux density 
of 11.138 Jy. We note that the final images were made with the 
entire data set. 

The c alibration of the data were carried out using standard 
Miriad dSault etalJI 19951) tasks. ATL0D was used to convert the 
raw data from RPFITS format to the native Miriad format. For 
the further calibration, both channels at the end of the bands 
were discarded due to a significant lack of sensitivity, hence a 
data set containing two times 13 channels with 8 MHz band- 
width each was obtained, yielding a net bandwidth of 208 MHz. 
Fortunately, our data do not suffer from self-interference typi- 
cal for the ATCA. These so-called "birdies" occur within the 
samplers at integer multiples of 128 MHz. The 19'' ! harmonic at 
2.432 GHz fortunately lies exactly between our two bands and 



Table 1. Parameters of the observations. Given are the number 
of pointings per field, the average integration time per pointing, 
the total area covered by the pointings, the rms noise and the 
resolution of the final images. 



Field 


N pn , 


tint 


area 


rms noise 


resolution 






h 


deg 2 


pjy/beam 


arcsec 


ELAIS-S1 


76 


6.17 


2.77 


70 


33.56 x 19.90 


CDF-S 


88 


4.13 


3.63 


80 


57.15 x 22.68 



Notes. For comparison, we note that the corresponding typical rms 
noise levels of the 1 .4 GHz images are 3 1 /jJy and for 38 yu Jy the ELAIS- 
Sl and CDF-S, respectively. 



therefore does not affect our dat a. After bandpass- calibration, 
RFI was removed using Pief lag (Middelbere 2006), a tool that 
compares statistics from channels with no or little RFI contami- 
nation to the other channels in the data and looks for outliers and 
sections with high noise. Bandpass calibration was carried out 
using the observations of PKS 1934-638. Phase calibration was 
done using the regular scans of the secondary calibrators, and 
then amplitude calibration was carried out using PKS 1934-638. 
The data were split into the individual pointings which were im- 
aged separately using uniform weighting. This provides a higher 
resolution at the cost of lower sensitivity than natural weight- 
ing. The cell size was set to 6" (3"). Because of the low frac- 
tional bandwidth (< 10%), a normal CLEAN procedure was used 
(in contrast to the Miriad multi-frequency clean implementa- 
tion MFCLEAN for the 1.4 GHz data) for deconvolution. Also the 
dynamic range in the pointings was sufficiently low as to not 
require any special deconvolution procedures. Cleaning was car- 
ried out using 2000 iterations, after which only marginal side- 
lobes around the two strongest sources remained, and the ef- 
fects of clean bias should be negligible. After cleaning, the mod- 
els were convolved with a restoring beam of 33.56" x 19.90" 
(57.15" x 22.68") at a position angle of -1.32 deg (-1.90 deg). 
The final mosaics were produced using Miriad's LINMOS task to 
account for primary beam attenuation and overlapping regions of 
the individual pointings. As a result, the edges of the final images 
exhibit a significantly higher noise level than the central regions 
(see Fig.Q]). 

3. The source catalogue 

Source extraction from the final 2.3 GHz images was performed 
in two steps: First, all known 1.4 GHz sources from ATLAS 
DR2 (Hales et al. in prep.) were used as an input catalogue to 
search for 2.3 GHz counterparts. This catalog contains a total of 
2240 individual 1.4 GHz sources across both the ELAIS-S1 and 
the CDF-S of which 631 could be reliably detected at 2.3 GHz. 
To identify these counterparts, the 2.3 GHz image was inspected 
by eye at the location of every individual 1 .4 GHz source. This 
was done by half-automatically browsing through an input cat- 
alogue, simultaneously displaying the corresponding regions in 
the 1.4 GHz and 2.3 GHz images and, if a 2.3 GHz counterpart 
was found, fitting a Gaussian to it and recording all fit parameters 
in an output catalogue. Thus, no blind search was conducted at 
2.3 GHz, allowing the inclusion of sources fainter than the typ- 
ical 5 <x cutoff in the catalogue. The lowest signal-to-noise ratio 
(SNR) among the 2.3 GHz sources is 3.9. 

In a second step, we checked for "2.3 GHz-only" sources. 
Those sources are unlikely to occur because of the noise levels 
of the 1.4 GHz and 2.3 GHz images. Since the faintest 1.4 GHz 
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Fig.l. 2.3 GHz noise maps of both ATLAS fields: CDF-S (left) and ELAIS-S1 (right). Contours start at TOyuJybeairT 1 (dark blue 
contour line) and increase by a constant value of 40yuJy beam -1 . 



sources in the DR 2 catalogue are as faint as 100 fjjy, whereas the 
faintest sources in 2.3 GHz are three times brighter, a 2.3 GHz- 
only source would have a highly inverted spectral index of a > 
2.4. Also, we note that, due to the temporal offsets between the 
1.4 GHz and 2.3 GHz observations, possible variability of some 
radio sources may lead to flawed spectral indices. Nevertheless, 
we used MIRIAD's SAD task to automatically produce a 2.3 GHz 
catalogue with a detection threshold of 5cr and matched the 
sources obtained that way to the 2.3 GHz catalogue obtained 
by cross-matching to 1.4 GHz sources. All sources identified by 
SAD could be matched to sources previously found as 1 .4 GHz 
counterparts, hence not a single 2.3 GHz-only source was found. 

To compute spectral indices, the flux densities in 1.4 GHz 
had to be re-measured because the resolutions of the 1.4 GHz 
and 2.3 GHz images differ significantly. Whilst the 1.4 GHz ob- 
servations included long baseline configurations of the ATCA 
(e.g. the 6A and 6B configurations), the 2.3 GHz observations 
were obtained with the ATCA in 750B and 750C configura- 
tions only. This leads to a resolution three to five times lower 
at 2.3 GHz compared to the ~ 10" resolution of the 1.4 GHz ob- 
servations, despite the shorter wavelength. Therefore, in order 
to avoid resolution effects, we convolved the original 1.4 GHz 
images with a Gaussian kernel to match the resolution of the 
2.3 GHz images as given in Table [2] We then used the ATLAS 
DR2 catalogue again as input to re-measure the 1.4 GHz flux 
densities in the low-resolution images. Because the image sensi- 
tivity was reduced by a factor of ~ 2 by this convolution (which 
is equal to tapering the data in the Mv-plane, so in fact weighting 
down long baselines), all sources in the low-resolution image 
should be a true subset of the hi-res sources and, therefore, all 
low-resolution sources should be detected when using the high- 
resolution catalog as input sample. 

During the re-measurement of the 1 .4 GHz flux densities us- 
ing convolved images, a problem surfaced with the ELAIS-S1 
image. Unlike with the CDF-S convolved 1.4 GHz image, where 
the procedure outlined above went well, the MIRIAD task imf it 
overestimated the flux densities by a significant amount (up to 
50%). This became obvious by inspecting residual images for 
sources with various flux densities and degrees of compactness. 
While the residuals looked good in the CDF-S, they showed 
bowls of negative flux in the ELAIS-S1 image. The cause of 



the effect rem ains unclear. The fitting routines in the software 
package CASA ([Jaegerl200 8) yielded similarly overestimated flux 
densities, so it appears that the effect is not caused by the fit- 
ting software, but is inherent in the images, and we had to work 
around this issue. This was done by iteratively fitting Gaussians 
to every individual source. The first Gaussian fit (which was al- 
ways overestimating the flux density) was used as prior. The re- 
sulting residual image was then examined in terms of total flux 
and rms noise. If the total flux in the residual image (trimmed 
to a size twice as large as the restoring beam) was negative and 
the rms noise was deviating from the median rms in the region 
of the source by more than 10%, the fit was discarded and re- 
placed by a fit with a peak flux density 2% lower than in the 
previous fit. The residual image was again examined in a simi- 
lar way. This procedure was iteratively repeated until the criteria 
(non-negative total flux and rms noise consistent with the rms in 
this region) were entirely met. As a last step, all final residual 
images were inspected by eye to assure the quality of the fit. 

A further problem in calculating spectral indices for ev- 
ery 1 .4 GHz was introduced by the very different resolutions of 
the 1.4 GHz and 2.3 GHz observations. Because the resolutions 
differ by a large factor, as explained above, numerous sources 
which are clearly separated in the high-res images are blended 
together in the low-res images. Therefore, we chose to com- 
pletely exclude such blended sources from our final catalogue 
if they could not be separated by simultaneously fitting two or 
more Gaussians. We note that this may bias our sample towards 
fainter sources since classic radio doubles (which are mainly ex- 
cluded by the blending) are generally brighter than single, iso- 
lated sources. 

The final catalogue was then constructed by merging the in- 
dividual catalogues from ELAIS-S1 and CDF-S. This resulted 
in a final catalogue comprising 63 1 radio sources with clear de- 
tections and well-measured flux densities in both 1.4 GHz and 
2.3 GHz. A portion of the catalogue is shown in Table [2] for il- 
lustration. 
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Table 2. A representative portion of the ATLAS 2.3 GHz source catalogue of the ELAIS-S1 field for illustration purposes. The CDF-S catalogue contains the same information. 



ID" 


RA (J2000) 


DEC (J2000) 


^ 1 .4 (jriz.ni 


Ail | A f • 1 1 . . 


.1 1 A I'l l . 1„ 


tW I A I'l l . 1„ 

1 -iu 1 .4 Lirlz.lo 


* 1 A UriZ.lO 


S2 H 


'-iu 2.3 Liriz 




2 1 

"1.4 


Aff" 

u l.4 


3 


class^ 




deg 


deg 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 










254 


51.770792 


-28.650297 


1.260 


0.084 


1.230 


0.093 


0.069 


0.561 


0.080 


0.075 


-1.58 


0.11 






178 


51.621458 


-27.427269 


2.740 


0.145 


2.220 


0.165 


0.122 


1.519 


0.107 


0.075 


-0.76 


0.11 






165 


53.319708 


-27.941922 


1.131 


0.054 


0.762 


0.158 


0.153 


0.828 


0.084 


0.073 


0.17 


0.13 


0.685 




163 


52.470667 


-27.492661 


2.160 


0.116 


2.070 


0.985 


0.980 


1.810 


0.116 


0.073 


-0.27 


0.95 






102 


52.561167 


-28.566331 


0.194 


0.037 


0.162 


0.155 


0.154 


0.344 


0.074 


0.072 


1.52 


0.64 






290 


53.190750 


-28.520756 


0.714 


0.044 


0.598 


0.126 


0.123 


0.400 


0.075 


0.072 


-0.81 


0.14 


0.211 




201 


52.399667 


-27.950881 


1.420 


0.087 


1.550 


0.153 


0.132 


1.154 


0.093 


0.072 


-0.59 


0.06 


0.067 


SF 


150 


52.122000 


-28.030825 


1.810 


0.087 


1.770 


0.281 


0.267 


1.016 


0.088 


0.072 


-1.12 


0.30 


0.902 




264 


52.105750 


-27.745814 


0.920 


0.054 


0.745 


0.176 


0.173 


0.674 


0.080 


0.072 


-0.20 


0.46 






323 


52.808292 


-28.785917 


0.834 


0.057 


0.818 


0.160 


0.155 


0.544 


0.077 


0.072 


-0.82 


0.38 


0.140 


AGNe 



Notes. The columns are as follows: (1) unique ID (2) degrees of J2000 right ascension (3) degrees of J2000 declination (4,5) 1.4GHz integrated flux density and flux density error as given in Hales 
et al. (in prep) (6,7,8) integrated flux density, flux density error as calculated using Eqn.[TJand local rms noise as measured from the low-resolution 1.4GHz image (9,10,11) integrated flux density, 
flux density error as calculated using Eqn.Q]and local rms noise as measured from the 2.3 GHz image (12,13) s pectral index calcu lated from the low-resolution 1.4GHz flux density and the 2.3 GHz 
flux density and its error as calculated by Eqn.|4](14, 15) spectroscopic redshift and spectral classification from lMao et alj d2012h if available. The entire table is available via the on-line version of 
the journal. 

<D) Unique only for the respective field (ELAIS-S1 or CDF-S). ^ Spectroscopic classification from M ao et al ]dioii). SF indicates an optical spectrum dominated by star formation activity while 
AGN indicates an optical spectrum typical for AGN (e.g. AGNe for an AGN with characteristic emission lines). E indicates that the spectrum is that of an early-type galaxy and hence the radio 
emission must originate from an AGN hosted by it. 
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We give an identifier for each source, the position, the orig- 
inal 1 .4 GHz flux density as measured by Hales et al. (in prep.), 
the re-measured 1.4 GHz flux density using the convolved im- 
age, the 2.3 GHz flux density as measured during the cross- 
matching process and finally the spectral index calculated from 
the matched-resolution 1.4 GHz and 2.3 GHz flux densities. All 
flux densities are integrated quantities. 

Errors for the flux densities were either taken from Hales 
et al. (in prep.) in c ase of the original 1 .4 GHz flux densities or 
calculated following lHopkins et al.l(l2003l) for the low -resolution 
1 .4 GHz flux densities and the 2.3 GHz flux densities, taking into 
account both the local rms noise at the source position as well as 
the quality of the Gaussian fit: 



| 2 2 
AS J L 1 image I 1 fit 

where [image quantifies the uncertainty in flux density due to 
the image noise and /i /,> represents uncertainties due to the fitting 
of a Gaussian. These tw o quantities can be calculated following 
IWindhorst etafld 19841) . 



ft image 



+ C 2 , + C 2 (2) 



2 ~/ ~P 

peak 



where <x is the local rms noise, Cf the error in absolute flux 
calibration and C p the error in flux due to eventual pointing er- 
rors. We here adopt a conservative val ue of C 2 + C 2 = 0.05 2 . 
The fitting error is calculated following ICondori(fl997l) . 



* + 1 



peak 



/ &B®b \ 


(A 


u 2 \ 

t*m 


\ ©M®m / 




©2 



(3) 



Here, &B®b are the restoring beam's major and minor axes, 
whereas ©m®»i represent the major and minor axes of the fitted 
Gaussian. The errors hs,Hm and fi m stand for the errors of the 
fitted Gaussian's peak flux, major and minor axes. The error of 
the spectral index was finally calculated using simple Gaussian 
error propagation, 



a -i II AS 1.4GHz \ 2 /AS 2.3GHz \ 2 
\\ J 1.4GHz / \ o 2.3GHz / 

where the factor of 2 « |[log(l. 4 GHz/2.3 GHz) ■ ln(lO)]- 1 ! 
arises from deriving the spectral index as defined in section Q] 
with respect to the flux densities, which are the only quantities 
with uncertainties in this equation (the error in frequency is ne- 
glected). 

Histograms illustrating the distribution of flux densities in 
the final catalogues are shown in Fig. [2] the relative flux den- 
sity errors as a function of flux density are presented in Fig. [3] 
Further spectral index properties are investigated in the next 
Section. 
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Fig. 2. Dis t ributio n of 1.4 GHz (hi-res) and 2.3 GHz flux densities in the final catalogue. AGN classifications of the sources are from 
Mao et al. I d2012h . 



4. Spectral index properties 

The distribution of spectral indices for our sample of 631 
1.4 GHz sources with unambiguously identified 2.3 GHz coun- 
terpart is shown in Fig. |4] 

With a median spectral index of -0.74 our sample matches 
the canonical value for synchrotron radiation of -0.7 very well. 
However, we note that the distribution is broad with a median 
absolute deviation of 0.35. We here use the median absolute de- 
viation defined as MAD = medianjOxj - median^ )|) as a mea- 
sure for the spread of a distribution since it is more robust against 
outliers than, for instance, the normal standard deviation, cr. In 
case of a perfect Gaussian distribution, cr and MAD are related 
via cr =s 1 .4826 MAD. This als o compares nicely to other studies 
such as Prand oni et al.1 (120061) whos e sample has MAD = 0.33. 
The spread of the llbar et all (120091) sample is slightly smaller 
with MAD = 0.27. This difference is due to the larger number of 
sources with flat or inverted spectra in our sample. We repeated 
the analysis including upper l imits (following the surviv al anal- 
ysis techniques developed in Fei gelson & Nelsonlll985l) for all 
1.4 GHz sources which do not have a 2.3 GHz counterpart and 
are not affected by resolution effects (see Sect. [3}. These sources 
do have a 2.3 GHz flux density less than three times the rms noise 
at the respective position. The median spectral index then be- 
comes slightly steeper with -0.81 and a MAD of 0.38. The main 
reason for only having this slight change is the large difference 
between the 1.4 GHz sensitivity (about 30^Jy) and the 2.3 GHz 
sensitivity (70-80/Jy). The number of sources detected at both 
1.4 GHz and 2.3 GHz is 631 whereas there are 552 sources with 
an unambiguous upper limit at 2.3 GHz. We note that the errors 
of our spectral indices become quite large for low flux densi- 
ties. This is because currently available radio data and the flux 
densities extracted from them have typical uncertainties at the 
10% level. This is mostly due to either inhomogeneous data pro- 
cessing and calibration strategies as well as the widely adopted 
method of fit ting 2-dimensi onal Gaussians to extract flux densi- 
ties (see e.g. ICondonlll997l) . Future large-scale continuum sur- 
veys, such as ASKAP/EMU, are attempting to reach flux den- 
sity errors of only 1 % by developing both novel, full y automated 
calibr ation as well as source extraction techniques dNorris et al.l 
1201 lbl) . 




io 1 io° io 1 io 2 



flux density [mjy] 

Fig. 3. The fractional flux density errors for all ELAIS-S1 and 
CDF-S sources as function of their flux density. The solid red 
line indictates how the fractional error should decrease with in- 
creasing flux density in an image with an overall rms noise of 
70 //Jy according to Eqn.|2] A calibration error of 5% is assumed. 



4. 1 . Spectral index vs. flux density 

A plot of the spectral index of our sources against their flux 
density is shown in Fig. [5] We point out that the sensitivities 
of our 1.4 GHz and 2.3 GHz data differ by a factor of about 
3. Therefore, progressively negative spectral indices cannot be 
measured anymore towards fainter 1.4 GHz fluxes (as indicated 
by the black dashed line in Fig. |5}- For instance, a source with 
S 1.4GHz = 200 ^Jy can only have a 2.3 GHz counterpart when 
exhibiting a spectral index ^ 0.5 since the 3<x detection level of 
our 2.3 GHz data is about 250^/Jy (see Sect. |2j. To account for 
this fact, we also included upper limits on spectral index for all 
1.4 GHz sources that unambiguously show no 2.3 GHz counter- 
part (the upper limits were calculated using $2.3GHz < 3c where 
cr is the local rms noise for a point source) and are not affected 
by resolution issues (see Sect. [3] e.g. a close double in 1.4 GHz 
that would be inseparable at the 2.3 GHz resolution). With these 
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Fig.4. Distribution of 1.4GHz to 2.3GHz spectral indices in 
the final catalogue for all sources with unambiguously detected 
2.3 GHz counterpart and both the 1 .4 GHz and 2.3 GHz flux den- 
sities exceeding 0.3 mJ y. AGN classifications of the sources are 
from lMao etafl (l2012h . 



upper limits, we calculated median spectral indices for differ- 
ent flux density bins to investigate a potential flattening of the 
mean spectral index towards fainter flux densities. As expected, 
the median spectral indices for only the detected sample (so the 
631 sources with 2.3 GHz counterparts) show a clear trend of 
flatter spectral indices with lower flux density (Fig. [3] blue cir- 
cles). Including the upper limits into the median calculation, this 
trend vanishes immediately (Fig. [5] blue diamonds): the lowest 
flux density bin (about 300-600 yidy) is completely dominated by 
the upper limits and can hence not be included in further anal- 
ysis whereas the two flux density bins at 700 fiiy and 1 mJy are 
now a even a little steeper than the ones at higher flux densi- 
ties. Hence, we find no evidence for a flattening of the spectral 
index towards the lowest flux dens ities. This is consistent with 
the findings bv lRandall et alJ (1201 2l) who used three frequencies 
(840 MHz, 1 .4 GHz and 2.3 GHz) for their spectral in dex calcu- 
lation and extend their findings to fainter flux limits. ITbar et al] 
(2009) come to the same conclusion, too. 

A flattening of t he spectral index wi th lower flux density lev- 
els, as_e ; g^_se^n_by|^ang|etal] ([2003), was initially suspected 
by IWindhorst et all il993 . A major line of argument is that, 
at sub-mJy levels, a change in population takes place: At mJy 
levels, the radio sky is completely dominated by AGN (see e.g. 
iPadovani et al.ll201ll) . so the dominant emission mechanism is 
non-thermal synchrotron radiation. An intrinsically flatter spec- 
trum would e.g. be produced by dominating emission due to 
thermal bremsstrahlung, but this ma kes up only a marginal frac- 
tion even in starburst galaxies (see ICondonlll992l) . In addition, 
recent work by several authors indicate that the transition from 
AGN-dominated to SF-dominated radio sourc es takes place 
below a flux density of 100/Jy. For instance, Padovani et al. 
d201 ll) argue that A GN still dominate th e source population 
down to 100/iJy and ISmolcic et al] d2008l) claim that the frac- 
tion of SF-powered radio sources has a nearly constant value of 
about 40% for flux l evels between 50 ^Jy a nd 700 /zJy. A similar 
result is obtained by Seymou r et al.l 0008) investigating the dif- 
ferential number counts in the 13 Hour XMM-Newton/Chandra 
Deep Field Survey. Hence, this change in the faint radio popula- 
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Fig. 5. Spectral index vs. flux density. The gray 2D histogram 
shows our sample of 631 sources with 1.4 GHz and 2.3 GHz de- 
tections. Since the sensitivity of the 2.3 GHz observations is ~ 3 
times lower than that of the 1 .4 GHz data, the black dashed line 
indicates the minimum spectral index that could be measured 
with our data. Therefore, the median spectral indices (blue cir- 
cles, x-error bars indicate the flux density bin size, y-error bars 
correspond to the MAD) for just the 2.3 GHz detected sample 
show significantly flatter values towards lower flux densities. 
When including also upper limits from the 2.3 GHz non-detected 
sources, this trend completely vanishes (blue diamonds). Hence 
we conclude that there is no evidence for a flattening of the mean 
spectral index with lower flux density. 



tion cannot account for our findings since we do not reach such 
faint flux density levels. 

With our catalogue of 63 1 sources with spectral indices, we 
now focus on the application of the spectral index in terms of 
its predictive quality. Since it is a quantity which is relatively 
easy to measure, the spectral index would be a useful tool to 
characterise sources in large-scale surveys, if a dependence on 
other properties of interest such as source type, redshift or other 
scaling relations can be established. 

4.2. Spectral indices as object type discriminator 

The radio spectral index is commonly regarded as a possible 
tool for discriminating differen t radio source types, hence emis- 
sion mechanisms. For instance. |Huvnh et al .1 d2007l) select AGN 
based on the assumption that, at low frequencies, only AGN 
can produce flat or inverted spectral indices because of self- 
absorption in an optically thick medium. Hence they propose a 
cut between a = -0.3 and a — to classify sources as AGN. 

With our source sample, we can examine this classification 
scheme. We cross-matched all our sources with and without un- 
ambiguous 2.3 GH z counterpa r t to th e spectral classifications for 
ATLAS sources bv lMao etalJ(l2012h . This cross-match lead to 
132 resultant sources of which we can identify 22 SF galaxies 
and 98 AGN among our sample with measured spectral indices. 
The remainder 12 sources could either not be reliably classified 
by means of their optical spectrum (4 sources), or they only have 
spectral index upper limits (5 AGN, 3 SF galaxies). We suspect 
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Fig. 6. Histogram illustrating the distribution of spectral indices 
of spectroscopically classified AGN and SF galaxies. 



that the reason for most spectroscopically classified sources hav- 
ing actual spectral index measurements is due to the fact that the 
target sources for spectroscopy were chosen such that stronger 
radio sources were favored against weaker ones. Therefore, the 
spectroscopically classified sample shows a significantly higher 
mean flux density at 1 .4 GHz so that most of these sources are 
detected at 2.3 GHz, too. Hence, since the vast majority of cross- 
identifications both have reliable spectral classifications and ac- 
curately measured spectral indices, we chose to omit the few 
sources that miss some information and proceed with these 120 
sources. A histogram showing the distribution of spectral indices 
among those 120 sources is presented in Fig. [6] 

As one can clearly see, the AGN population shows signifi- 
cantly flatter spectral indices than the SF galaxies with median 
values of -0.63 for AGN and -0.85 for SF galaxies. Note that 
there are only two SF galaxy exhibiting a spectral index greater 
than -0.5. Therefore, fiat or inverted spectral indices can reliably 
be associated with dominant AGN activity. On the other hand, 
steep spectral indices (i.e. a < -1.0) are not associated with any 
source type as the histogram clearly shows. Therefore, uniquely 
steep spectra cannot be used as discriminator for source type. 

Because the vast majority of (radio) sources detected in con- 
temporary and forthcoming surveys will not have any spectro- 
scopic information, we also searched for correlations between 
the spectral index and other photometric object type indica- 
tors. A popular way of photometrically selecting AGN with- 
out being biased by extinction is the use of an infrared color- 
color diagram based on the four channels of the IRAC instru- 
ment (Fa zio etaD 120041) aboard th e Spitz er Space Telescope. 
Originally proposed by lLacv et all (120041) and therefore often 
referred to as the "Lacy diagram", this plot has bee n further 
investigated (for an overview see ISoifer et all |2008) so that 
there are currently reliable template spectral energy distribu- 
tions (SE Ds) available covering the entire optical to MIR regime 
(e.g. lAssef et~aDl2010h . To construct a "Lacy diagram" for our 
source sample, we cross-matched the radio sources with the 
SWIRE catalogue, data release 3, to obtain IRAC flux densities. 
We found 265 secure identifications where all four IRAC fluxes 
are available. The corresponding plot is presented in Fig. [7] For 
comparison, we also plotted all SWIRE-detected sources in the 
CDF-S and the ELAIS-S1 field as color-coded, 2-dimensional 
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Fig. 7. "Lacy diagram" of our source sample. For comparison, 
all sources detected in all four IRAC bands within the SWIRE 
survey of ELAIS-S1 and CDF-S are plotted as 2-dimensional 
histogram. The black dashe d line indicates the l ocus of AGN in 
this diagram as proposed bv lDasvra et al.l d2009h . 



histogram. The locus of potential AGN in this diagram with re- 
spect to various contaminating effects, predominantly line emis- 
sion by p olycyclic aromatic H ydrocarbons (PAHs), was deter- 
mined by Das vra et al.l (2009) and is indicated in Fig. UJ with a 
dashed line. We split our source sample in three sub-samples 
based on their spectral index: steep spectrum sources are classi- 
fied as having a < median(a) - MAD, normal spectrum sources 
satisfy median(a) - MAD < a < median(o') + MAD and 
flat/inverted spectrum sources show a > median(o') + MAD. 
One can see that flat or inverted spectrum sources are more fre- 
quent on the "AGN branch" of the diagram whereas the loca- 
tion of star-forming galaxies (the branch that goes straight up) 
is more or less avoided by these sources. This "star formation 
branch"(log(5 5. 3.6^111) < 0) is dominated by normal and 
steep spectrum objects with only a very few flat/inverted spec- 
trum sources. 

To summarise, the radio spectral index can at least be used 
to identify a fraction AGN, namely those with optically thick 
media that cause their radio emission to flatten because of self- 
absorption energy losses. However, spectral indices are not a 
strong discriminator for the source type because (i) outliers from 
this general trend are frequent and (ii) optically thin AGN can- 
not be separated from SF galaxies with spectral indices. This is 
further highlighted in Fig. [8] where we plotted the "Lacy dia- 
gram" for all A TLAS sources with spectroscopic classification 
dMao et al1l2012h . It becomes obvious that the "AGN branch" in 
this plot is really occupied by AGN only whereas the "star for- 
mation branch" both hosts SF galaxies and AGN. Hence, also 
a selection based on MIR color criteria leads to the result that 
a AGN can be reliably selected without contamination by star 
forming objects. The other way round, selecting SF galaxies 
only without AGN contamination, is not possible by means of 
MIR color selection, similar to a selection based on radio spec- 
tral indices. 
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Fig. 8. "Lacy diagr am" of our source sample with spectroscopic 
classification from Ma o et all (120121) . As in Fig. [7] all sources 
detected in all four IRAC bands within the SWIRE survey of 
ELAIS-S1 and CDF-S are plotted as 2-dimensional histogram. 
The black dash ed line indicates the locus of AGN in this diagram 
as proposed bv lDasvraet al.l (|2009). 



4.3. The z-a relation 

The z-a relation is an empirical relation seen between 
the spectral index of radio sources in a flux-limited sam- 
ple and their redshifts. It states that higher-redshift radio 
sources tend to exhibit steeper spectral indices, typically with 
a < -1. It was found by various authors using different 
large-scale radio s urvey s at the Jy and mJy level, see e.g 
Athr eva & Kapahil dl998l).|Pe Breuck et al] d2000tl200lL [2006) 



Klam er et al.l (l2006h . lAfonso et all (1201 ll) . and lKer et al.f(l2012l ) . 
Th e first quantitative formulation of this relation was given 
by IVerkho danov & Khabibul linal (1201 Ol) using a large sample 
(~ 2500) of distant radio galaxies selected from various sur- 
veys across the entire sky. Above a redshift of 1, they find clear 
evidence for such a correlation and quantify it by a linear fit 
a = (-0.70 + 0.02) + (-0.15 + 0.01)z out to z ~ 5. The currently 
favoured explanatio n for this effect is based on cosmic expansion 
dKlamer et al.l l2.006): At higher redshifts, when the universe was 
denser on average than at present, the radio jets of high-z AGN 
had to penetrate into a denser intergalactic medium where they 
are compressed, similar to the jets of radio galaxies located in 
the cores of dense clusters in the nearby universe. This compres- 
sion then leads to the characteristic steepening of the spectral 
index due to synchrotron and inverse-Compton energy losses of 
the relativistic electrons. 

With the spectroscopic information from iMao et al.l (1201 2l) . 
we have investigated the z-a relation in our sub-mJy sample. 
From our initial 631 sources with measured spectral index, 169 
(27%) have spectroscopic redshifts. In addition, we have 236 
sources with spectroscopic redshift and an upper limit on a only. 
A plot of these redshifts versus the spectral index or the up- 
per limits, respectively, is shown in Fig. [9] We again note that 
the 169 sources with spectroscopic redshifts and actual spec- 
tral index measurements are not affected by sensitivity bias (see 
Sect. 3]l since they all have 1 .4 GHz flux densities in excess of 
0.8 mJy and hence spectral indices of at least -2.0 can be mea- 
sured for them, given our 2.3 GHz sensitivity. 



Fig. 9. Redshift versus spectral index plot for all sources with 
spectroscopic redshift in our sample. Sources with an actual 
measurement of a are split in spectroscopic classes (AGN: black 
square, SF: black star, unknown: black diamond). Sources with 
an upper limit on a only (black dots) were not split up. Note that 
R 2 = 0.09 for the fit to our data points (black symbols) indicates 
highly uncorrelated quantities. A slight correlation is usually as- 
sumed when R 2 > 0.5 with good correlation for R 2 > 0.85. 



Surprisingly, in this diagram, there is no obvious correla- 
tion between spectral index and redshift. Note that this result 
is also independent of the spectroscopic classification of the ob- 
jects. But since radio sources with significant redshift are, due to 
the sensitivity limits of currently available telescopes, all AGN- 
powered (even a galaxy forming stars at a very high rate of 
1000 M yr~' would be undetectable within ATLAS for any red- 
shift z > 1 .2), the z-a relation is yet only investigated for AGN. 

This finding (the absence of a z-a relation) contradicts the 
contemporary literature, and we offer two potential explanations: 
First, our sample of sources is very small (only 9 sources) for 
the redsh ift range z > 1 were the z-a relati on becomes impor- 
tant (see Verkhoda nov & Khabibullinall2010l) . Hence, our find- 
ings rely on small number statistics and have therefore to be 
treated with care. Second, there could be a change in the radio 
source population between our sub-mJy sample and the sam- 
ples used in other work which are mainly much more luminous 
sources with flux densities of several tens to hundreds of mJy 
(despite their significant redshifts , so these sourc es are very lu- 
minous). We also point out that iKer et al] (1201 2l) found only a 
slight correlation between spectral index and redshift. They ar- 
gue that the relation becomes fully applicable only for resolved 
sources, if the explanation is correct that the z-a relation arises 
from high-redshift radio lobes working against the denser inter- 
galactic medium. This reasoning is supported by our findings 
because nearly all of our sources are compac t. In order to ef- 
fectively select potential high-z radio sources, iKeretal] d2012l) 
suggest to use the relation between Tif-band flux density (or also 
Spitzer 3.6 yt/m IRAC channel) and redshift which is much tighte r 
than the traditional z-a relation, see also Norri s etaLl (120 1 la). 
This relation was also investigated by Middelberg et al. (2011) 
to select a certain class of objects which is suspected to be lo- 
cated at very high redshifts (z > 2), the so-called Infrared-Faint 
Radio Sources (IFRS). 
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Fig. 10. S 20cm/S 3.6fim f ux ratios for our spectral index source 
sample. Note the clear but not very tight trend of increasing red- 
shift with increasing S20cmfS3.6fim- The limiting value for IFRS 
is shown as red dashed line. 



According to IZinn et al] d201 ll) . a source is deemed an 
IFRS if its measured ratio between radio (20 cm) and near- 
infrared (NIR, mostly Spitzer IRAC 3.6 fim) flux density exceeds 
a cutoff-value of 500 and, at the same time, its NIR flux den- 
sity is less than 30/iJy. Those extreme sources (often exhibit- 
ing no 3.6 yum counterpart at all, even in ultra-deep or stacked 
NIR images with noise levels reaching the lOOnJy level) are 
very difficult to investigate, in particular regarding their redshift, 
since optical spectroscopy can only be done for such IFRS that 
at least have moderately bright optical counterparts. Therefore, 
we only have very few sources in our sample with spectroscopic 
redshifts that meet at least the first IFRS criterion. A real IFRS 
meeting both criteria is not in our sample. A plot illustrating 
the S l 20cm/>S'3,6j"m A ux ratios of our source sample with respect 
to redshift is shown in Fig. [10] Contrary to the z-a relation in 
Fig. [9] there is a clear trend of increasing redshift with increas- 
ing S20 1:111/53.6,™ flux ratio. Although the scatter in the relation is 
fairly large, in particular for high flux ratios, this relation seems 
to be a much better estimator for the redshift of a source than the 
classical z -a relation. This is in very good agreement with the 
findings o f lKer et al.ld2012l) who also see only a marginal corre- 
lation of the spectral index with redshift but state that by far the 
most efficient selection criterion for high-z sources in complete 
samples is based on /T-band magnitude. 

We conclude that our data do not confirm the z-a relation, 
at least for faint radio sources. This may have several explana- 
tions: On the one hand, the popular belief that the z-a relation 
is caused by the denser intergalactic medium (IGM) at higher 
redshifts against which the radio j ets of powerful sou rces have 
to work and thereby lose energy ( Kla mer et alj|20 06). may not 
be applicable to faint sources since they are mostly compac t. For 
instance, the simulations used by iRaccanelli et al.l (1201 ll) pre- 
dict that compact Fanaroff-Riley type I (FR I) sources are more 
frequent at any redshift greater than 1 by several orders of mag- 
nitude compared to extended FR II sources. Therefore, the over- 
all radio population, at least at the sub-mJy level, s hould mainly 
consis t of such compact FR I sources for which the lKlamer et al.l 
(2006) explanation is less suitable. Furthermore, the z-a relation 
may only be an apparent correlation. Nearly all radio spectra of 
the highest-z (z > 4) radio sources known to date are signifi- 



cantly curv ed or even peak ed at observed-frame frequencies of 
<500MHz dKer et al.ll2012L , their Fig. 5), a behavior that is in- 
dicative for either young sources or sources w ith recent merger 
activity (e.g. lO'Deall 19981: iRandall et alJl201lh . Hence, spectral 
indices calculated between 1.4 GHz and another radio frequency 
commonly in use today, e.g. 610MHz or 5 GHz, measure the 
high-frequency flank of such a peaked spectrum which yields 
increasingly steep spectral indices the more the radio spectrum 
is curved or peaked. Therefore, sources selected by a steep spec- 
tral index are predominantly of these two types (either young or 
merger). Because those sources are more frequent at higher red- 
shifts when the universe was denser (more merger) and younger 
(more young sources), it only appears that a selection based on 
steep spectral indices yields high-redshift objects despite the real 
astrophysical phenomenon that is selected is completely differ- 
ent, namely young or merger sources. 

4.4. Spectral indices and the radio-IR correlation 

The radio-IR correlation is deemed to be one of the tightest 
correlations in astrophysics. It states that the radio emission of 
a galaxy is positively correlated with its mid- and far-infrared 
emission due to recent star-formation activity. The IR emis- 
sion in this scenario is produced directly by the star-formation 
processes whereas the radio emission comes from supernova 
remnants related to the star- formation event. Af t er its discov- 
ery in the late 1980s (e.g. Ide Jong et all 119851: ICondon et ail 
1991), it was subjected to intensive studies while the tec hnol- 
ogy of infrared telescopes advanced ( Appl eton et alj 2004) and 
is no wadays proven to hold over a wide range in both flux den- 



( Huvnh et al 



2010tlMaoetalJl2011l) 



. . nge 
sitv dBovle eTall 120071: iGarn & Alexander! 120091) and redshift 



We here use co- located Spitzer obser vations at 24jum from 
the SWIRE survey dLonsdale et al.l 120031) to investigate the de- 
pendence of the correlation parameter, ^24 = log^^m/^Ocm) 
on the radio spectral index. Matching our initial sample of 63 1 
radio sources to the SWIRE IR catalogue with a match radius 
of 5", 145 counterparts were identified. Because the radio-IR 
correlation should only hold for star-forming objects, we ex- 
cluded all known AG N among these so urces based on the spec- 
tral classifications by iMao et all d2012l) and other AGN in dica- 
tors like morphology wh ich was done in iNorris et alj j2006) and 
Middelberg et al. (2008). After this rejection procedure, we were 
left with a total of 86 sources with secure 24 yum detection and 
no signs of AGN activity. A plot showing $24/™ versus S 20 cm 
for these objects is presented in Figure fTTI 

The correlation in our source sample is , at l east fo r their 
slope, well described by the fit of lAppleton et alJ d2004l) which 
was obtained using SWIRE data but over many fields, not just 
ELAIS-S1 and CDF-S. The few sources that lie at the top-left 
corner of the diagram, exhibiting an excess of radio emission, are 
leftover AGN which were not removed from the sample during 
the AGN rejection step. 

Calculating the correlation parameter ^24 for our sample, 
we find a median value of 0.68 with a fairly large spread of 
MAD = 0.36. When plotting the individual values of ^24 against 
the spectral indices of the sample (Fig.fTZli. we see no clear cor- 
relation between the two quantities. 

Nevertheless, a general trend could be found that sources 
with a spectral index steeper than the median value show a much 
larger spread in <?24 values than sources with a spectral index 
flatter than the median. This is increasingly surprising with the 
results of Sect. 14.21 where we found that flat or inverted spec- 
tral indices are indicative for AGN-powered sources. Since AGN 
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Fig. 11. Diagram illustrating the correlation between the flux 
densities at 24yum and 20cm (1.4 GHz). The sources are divided 
in three spectral index bins according to a < median(o') - MAD, 
median^) - MAD < a < median(a) + MAD and a > 
median^) + MAD where MAD is the median absolute devia- 
tion of the spectral index of the entire 24 /mi-detected sample. 
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Fig. 12. <?24 values of our source sample plotted against spec- 
tral index. The solid red lines indicate the median of the two 
quantities while the dashed lines indicate the spread in terms of 
(median+MAD). 



should not obey the radio-IR correlation, this tightening of the 
scatter in ^24 with increasing spectral index is still to be inves- 
tigate d. One possible exp lanation could arise from the findings 
of iPadovani et alj d201 lb who investigated star forming galax- 
ies, radio-loud and radio quiet AGNs in the CDF-S. They found 
that regarding their luminosity function and its evolution, radio- 
quiet AGN are pretty much similar to star forming galaxies and 
not to their radio-loud siblings. Although we aimed to exclude 
AGN from the sample discussed in this section based on indi- 
cators such as radio morphology or spectroscopy, we predomi- 
nantly eliminated radio-loud AGN and several radio-quiet AGN 
may still contaminating the sample (as suggested by the few flat 
and inverted spectrum sources in Fig.[T2l. That these radio-quiet 
AGN do obey the radio-IR correlation is a natural consequence 



of the findings of Padov ani et al.l (1201 ll) and already discussed in 
their work. 

In summary, we conclude that the radio-IR correlation is not 
affected by the radio spectral indices of sources in a flux-limited 
sample. This is yet another proof for their remarkable universal- 
ity which has been shown for redshift, flux density levels and 
now also for radio spectral properties. 

5. Summary and outlook 

We have conducted deep (<x = 70 - 80 fiiy) radio observations of 
the ELAIS-S1 and CDF-S regions at 2.3 GHz using the ATCA 
to supplement the sources detected at 1.4 GHz within ATLAS 
with spectral index information. Our final catalogue with sources 
well detected at both frequencies and therefore with high quality 
spectral indices comprise a total of 63 1 sources and is available 
in the online version of this article. The median spectral index 
of our distribution is -0.73 with a median absolute deviation 
(MAD) of 0.35. We have investigated the properties of the pop- 
ulation with respect to spectral index and tested correlations be- 
tween the spectral index and various other quantities. Our main 
findings are as follows: 

1. We do not find evidence for a flattening of the spectral in- 
dex towards fainter (several 100 fiiy) flux densities. This is in 
good agree ment with the majority of previous wo rk by other 
authors (e.g. iRandall et alJl20rl llbar et alJ[2009h . Since we 
also find that star-forming galaxies tend to have steeper ra- 
dio spectra compared to AGN, we can support the current 
notion that star-forming galaxies become the dominant type 
of radio sources only at flux densities below lOOji/Jy since 
our sample does not reach such low flux densities. 

2. Considering the spectral index as discriminator for the ob- 
ject type, we find that flat- or inverted-spectrum sources are 
mostly AGN while SF galaxies show significantly steeper 
spectra with a median spectral index of -0.85 (AGN median 
spectral index is -0.68). This flattening is caused by self- 
absorption losses within optically thick AGN under phys- 
ical conditions which are not in general met in SF galax- 
ies. However, there is a significant number of outliers which 
affects the selection of AGN based on spectral indices. 
Furthermore, optically thin AGN exhibiting the canonical 
spectral index of -0.7 for synchrotron radiation can not be 
separated from SF galaxies. Therefore, we conclude that the 
spectral index is not a strong discriminator for the source 
type but can be helpful if there are no other, more sophisti- 
cated diagnostics. 

3. We investigate the z-a relation which states that high- 
redshift radio galaxies show steeper spectral indices. The re- 
lation cannot be confirmed with our data. We suspect that 
this is explained either by the fact that the currently f avored 
explanation for the z-a relation dKlamer et alJl2006h is not 
applicable to the sub-mJy population or because of the com- 
plete failure of the z-a relation which may only be an ap- 
parent correlation. Since most high-z radio galaxies show 
curved or even peaked radio spectra, indicative for young or 
merger sources, a selection based on steep spectral indices 
may predominantly select those sources. But because both 
young and merger sources are more frequent at high redshifts 
when the universe was denser and younger, it appears to se- 
lect high-z sources. 

4. We find no evidence for a dependence of the radio-IR corre- 
lation on radio spectral index. This therefore widens the uni- 
versality of this correlation which is already known to hold 
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over a wide range of flux densities and redshifts. However, 
we do see a trend for the scatter of <?24 which is larger 
for sources with steep spectral indices and becomes smaller 
while going to positive spectral indices. 

In conclusion, we want to point out that radio spectral in- 
dices can be useful to characterise a radio-selected sample, but 
since there are only rough correlations with other source proper- 
ties, the predictive power of spectral indices may have been over- 
rated in the past. This is unfortunate, in particular for future all- 
sky rad io surveys such as th e Evolutionary Map of the Universe 
(EMU. [Norris et al.ll201 lbl) to be carried out with the Australia 
Square Kilometre Array Pathfinder (ASKAP) telescopes in the 
southern sky or similar surveys in the northern sk y planned to be 
done with the Low Frequency Array (LOFAR, Morganti et al. 

or the upgraded Westerbork Synthesis Radio Telescope 
when equipped with th e new focal plane array system Apertif 
(iRottgering et al.ll201 ll) . A significant fraction of the sources de- 
tected by these surveys at a typical noise level of 10 fjjy will 
have no other than radio data available, even c onsidering fu- 
ture surveys such a s Pan-Starrs dKaiser et alJl20ld) . SkyMapper 
( Kell er et alJ [2007). and the Large Synoptic Survey Telescope 
dlvezic et alj I200ST) in th e optical or the W ide-field Infrared 
Survey Explorer (WISE, IWright et all l2010h and ESQ public 
surve ys with the VISTA and VST telescopes fArnabold i et alj 
120071) in the infrared. 
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